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Abstract. In this article we prove local well-posedncss of quasilincar disper- 
sive systems of PDE generalizing KdV. These results adapt the ideas of Kcnig- 
Ponce-Vega from the Quasi-Linear Schrodinger equations to the third order 
dispersive problems. The main ingredient of the proof is a local smoothing es- 
timate for a general linear problem that allows us to proceed via the artificial 
viscosity method. 

1. Introduction. In this paper we consider the following system of PDE: 

d t u + a(x, t, u, d x u, d^u) ■ d^u + b(x, t, u, d x il, d^u) ■ d 2 x u 
+c(x, t, u, d x u) ■ d x u + d(x, t,u) ■ u = f(x, t) on Rx/ (1-1) 
u(x, 0) = uq(x) on R x{0} 

where / c [—1,1] is a time interval containing 0; all functions are real valued; 
u = (tii, • ■ • , u n ) is an unknown; u Q = (u ,i, . . . , «o,n) an d / = (/l, • • • , fn) are given 
data; and coefficients a, b, c, d are n x n matrix valued functions. 

This way (1.1) is a generalized KdV with the dispersive relation a that depends 
on the unknown and its derivatives (u, d x u, d x u). 

Moreover, well-posedncss of (1.1) allows one to study fully non-linear dispersive 
systems of the form: 

d t v + f(x, t,v 7 ... dlv) = for (x, ()elx [-T, T] 
v(x,0) = v (x) 

if the non-linear function / satisfies appropriate assumptions similar to the one for 
(1.1). Indeed, differentiating this equation and letting u = (v,d x v) allows one to 
reduce the problem to a quasi-linear one. 
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Well posedness of the semi- linear analogues of (1.1), where the top order a = 1, 
particularly when 6 = and c is a polynomial, is quite well-understood with Local 
Smoothing and Strichartz estimates playing a significant role, cf.. [13] and [15] and 
references within. While quasi-linear dispersive equations are of interest in physical 
applications, in particular to water waves with variable dispersion, they are far less 
understood. 

Well posedness has been established for quasilinear dispersive equations with 
special algebraic structure for which conservation laws have been found, for example 
the following shallow water wave equation 

Uf ^txx ~t~ SuU x — ^UxlLxx -\- UUxxx 

see [3] and references therein. However, finding such conservation laws is not pos- 
sible in general for (1.1). 

A major advance for the well-posedness of a scalar (1.1) was a work of Craig- 
Kappeler-Strauss in [5]. However, in addition to being restricted to scalar equation, 
they had to make a technical assumption for the favorable sign of the coefficient 
b, which was not natural in light of the semi-linear results in [15]. However, the 
pioneering work of Kenig-Ponce-Vega [14] for the Quasilinear Schrodinger equa- 
tion suggested a method to prove the well-posedness of (1.1) under more general 
assumptions on the coefficients. 

The method of Kenig et al. was a modification of the energy method, which was 
a successful approach to treat well-posedness of quasi-linear wave equation in high 
regularity Sobolev spaces in the 70's cf... [10]. Namely, the energy method relies 
on estimates of the form || u\\h° = O(||ito||ff s X which are proved via integration 
by parts, symmetry of the top order terms, Sobolev embedding and Grownwall 
inequality. However, the energy method cannot in general work for the (1.1) without 
modification due to the infinite speed of propagation. Overcoming this difficulty 
of controlling the size of a solution by the data occupies most of this paper. The 
heart of the matter can already be seen, when trying to prove L 2 well-posedness for 
a linear case of (1.1), regularized by a parabolic term for < e < 1. 

{d t u + a(x, t)dxil + b(x, t)d 2 u + c(x, t)d x u + d(x, t)u = —sdxU + f(x, t) ^ 
u(x, 0) = u (x) 

The standard energy method is available to prove an L 2 estimate for t > 0, only 
when b > 0. However, modifications of the energy method are needed, in general. 
Below we review the several works that motivated the approach used in this paper. 

In the case of a = Id, c symmetric and an integrable b, Kcnig-Staffilani in [15], 
motivated by [9], were able to cancel bd 2 term with a change of variables v(x,t) = 
&(x, t)u(x, t) for an appropriate After this change of variables a standard energy 
argument (as explained above) works for v and hence gives the L 2 estimate for 
I = [0, T] small enough 

\\u\\ L?Ll < A(||t?o(aO|| ia + \\f{x,t)\\ L \Ll) (1-3) 

where ||/(x, ^Hli^ = J r \\ f (t)\\ L2dt and likewise for all other space-time norms from 
now on. Similar argument was at the heart of the energy estimate that Lim-Ponce 
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used to prove well-posedness of a quasilinear Schrodinger equation in 1 dimension 
in [17]. However, this argument does not seem to work for a general system (1.2) 
with merely a symmetric top order ad^, or for the Schrodinger equation in more 
than one space dimension. As we show in section 2 the symmetry of the coefficient 
a is sharp for the estimate (1.3), in the sense that without this assumption, this 
estimate can be false. This suggested to proceed via a Local Smoothing estimate 
argument, with the side benefit of capturing the rcgularization effect for equations 
(1.2) and (1.1). 

Kenig- Ponce- Vega in [14] proved wcll-poscdncss of quasi-linear Schrodinger equa- 
tion in higher space dimensions by proving a Local Smoothing estimate generalizing 
from the case of a time independent variable coefficient linear Schrodinger in [7] 
and [6]. Local Smoothing for (1.2) means that for 8 > h 

IIt^jWsII^^^II^WIU' + IK^ 4 ^}" 1 /^,*)!!^) (1.4) 

(x) ' 1 

where we use the notation (x) = (1 + \x\ 2 )^ and interpret (d x ) u(£) = (£) fi(£) as a 
Fourier multiplier and hence (1.4) means that the solution of (1.2) is 1 derivative 
more regular than uo and 2 than / at a cost of weights. Note, that this effect is 
local, as the (1.2) is time reversible and by (1.3) cannot gain smoothness globally. 

Local Smoothing was first proven for KdV in [11], [16] and for the linear Schrodinger 
dtU + iLu = / with L = A in [4], [20], [21], [12] where the gain is \ derivative rel- 
ative to wo and 1 derivative relative to / respectively. This was generalized to 
C = a,ij(x)didj + 6(x)V in [6] and [7]. In [8] Doi showed that, roughly speaking, 
under appropriate asymptotic flatness of the coefficients for C above local smooth- 
ing is equivalent to the non-trapping of the bicharacteristic flow generated by the 
principal symbol of the operator C. Note that in the one spacial dimension, which 
is the relevant setting for this paper, the non-trapping condition is automatic by 
the coefficient assumption (NL1) and we will omit it. 

Our proof of the wcllposedness also involves using Local Smoothing to prove the 
energy estimate, however as (1.1) is of higher order than the Schrodinger equation 
the argument of [7] has to be modified. Using this method to prove wellposedncss 
of higher dispersive systems will be done in the subsequent work. 

To motivate the function space we use to prove well-posedness, we note that in 
order to prove the main linear estimate (1.3) a decay of the coefficient b is necessary. 
This phenomenon is similar to the Mizohata condition, which shows the necessity 
of sup X)t i w i =1 f Q ^sb(x + suj) ■ uods < oo for the L 2 well-posedness of Schrodinger 
equation dtu + iAu + b(x)Vu = 0, cf... [19], and we prove a corresponding result for 
(1.2) explicitly in the section 2. On the non-linear level of (1.1) this suggests an L 1 
condition on u and weighted Sobolev spaces provide a natural way to ensure an L 1 
condition in an L 2 -based Sobolev space. This motivates working with the weighted 
Sobolev spaces H s ' 2 for s G Z + , which we define as follows: 

k 2 s+3j 

j=0 j=0 a=0 
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Where \\f(x)\\ H » = ||(l + |f| a )*/(OIU' «. \\d s J(x)\\ L 2 + \\f(x)\\ L *. The space 
can be shown equivalent to H H s ((x) dx), which is a typical choice to ensure 

extra decay of the lower order term 69 2 as explained above. 

Coefficient assumptions. We state the precise assumptions on the coefficients a 
- d of the equation (1.1): 

(NL1) Dispersive property. The top coefficient a is symmetric, that is 

Oij(x, t, z) = dji(x, t, z) for all (x, t, z) € Dm eRx[-1,1] x [— M, M] 3n . 

Moreover, it is uniformly positive definite in Dm at time 0, that is for 
every M > there exists a constant Am > 

aij(x,Q,z)€i£j > A A/ |£| 2 

whenever (x, z ) e M x [-M, M] 3n . 
(NL2) Regularity. Let J € Z + be a given positive integer. We assume that all the 
coefficients a(x,t, z)-d(x,t, z) £ C\B J x Cl , with the space B J of C J functions 
with all J derivatives bounded. Specifically, there exists a function (J, M) — > 
Cj.m increasing in J and M > separately, such that for each J and M the 
coefficients 

a,,, fey e C}Bi s {D M ), en E C^^lxhM] x [— M, M] 2 "), 
with norms bounded by Cj.m, e-g- for a that means 

sup ||9 t Q 9f5>(x,i,z)|| L ~ <Cj.a/ 

0<a<l:0</3+|7|<J A/ 

(NL3) Asymptotic flatness and decay of linear parts. There exists 6' > \ and a 
constant Co, such that 

\lx) 2& ' d x a{x, 0, 0)|| L oo + \\(x) 2S ' b(x, 0, 0)\\ L oo < C 
\\(x) 25 ' d t d x a(x,t,0)\\ L ~ x ,+ Wix) 25 ' d t b(x,t,0)\\LZ? XI <C 
\\(x) S ' d x b\x,0,0)\\L°° + ||(.T) 5 'c t (a;,0,0)||L- < C 
\\(x) s ' d t d x tf {x, t, 0) Hi- 4- |Kx) 5 ' fleet (S) tj o) || irx; < Co 

where fet anc j c t are ^j^g antisymmetric parts of 6 and c respectively, defined 

While any 5' > | works in the (NL) assumptions, provided the definition of i/ 5,2 is 
modified with a weight (x) replaced by (x) S , for simplicity we set 8' = 1. 

Note, that unlike the constants Am in (NL1), which depend only on data uq, 
Cj t M in (NL2) is a family of constants that depend on the smoothness of the 
coefficients J and the size of the solution M. This is a natural assumption as we 
consider nonlinear equations and the coefficients may grow with the size of the 
solution. The size of the solutions is a priori unknown and is one of the quantities 
to be estimated. Thus when using (NL2) we will be explicit in the choice of J and 
M to avoid a possible circularity. 
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The proof is based on the artificial viscosity regularization of the initial value 
problem (1.1) for a parameter < e < 1: 

id t u" = -ed^-N s ^ + f 
\u £ (x, 0) = u (x) 

Where N u is the operator from (1.1) 

Ng = a(x, t, u, d x u, d 2 u) ■ + b(x, t, u, d x u, d x u) ■ d 2 

+ c(x, t, u, d x u) ■ d x + d{x, t,u) ■ I 

We then construct solutions u e for the regularized problem and show that the 
solution u — lim £ _j.o u E in the desired topology. 

Main Theorem. 

Theorem 1.1. Suppose the coefficients of (1.1) satisfy assumptions (NLl)-(NLS) 
with bounds Xm, Co and a function Cj^m for J > 19. Then (1.1) is locally well- 
posed in H ' . 

Specifically, by well-posedness we mean the following: 

(1) Existence. Let R > be given. Then there exists a T > 0, such that if the 
data (m , f) G Y = H s ' 2 x (l^_ hl] H s ' 2 n Cf^H^ , such that 

\\(u J)\\y = \\u \\ H s,2 + \\f\\ L i_ i i]H s.2 + ||/1| c p_^ i]i? 4,2 < R (1.8) 

then (1.1) has a solution u € X I: where / = [0, T] and Xi = C^H s ' 2 nCjH 4 ' 2 . 

(2) Uniqueness, u is the unique solution of (1.1) in X[. 

(3) Continuous dependence. The flow map (uq, f) — > u is continuous between Y 
and Xj. 

(4) Persistence of regularity. If (uq, /) & YH H s ' 2 x L^_ 1 ^H s ' 2 for s > 8, then u 
is in addition in CjH s ' 2 . 

We make a few remarks about the Theorem 1.1. 



We first note, that the time reversal u(x,t) — > u(—x,—t) applied to (1.1) pre- 
serves assumptions (NL1)-(NL3), thus Theorem 1.1 extends to [— T, 0]. The solu- 
tion we construct has enough regularity to be classical and the uniqueness proved 
above is valid for / = [— T, T}. Thus Theorem 1.1 produces a solution on [—T,T]. 
However, the parabolic regularization that we use to construct the solution requires 
a fixed sign of time and we stick to positive time intervals except for the uniqueness. 

The choice of well-posedness in H 8 ' 2 and the choice of T is determined by the 
main estimate Theorem 2.1 that proves (1.3) and its corollaries in the non-linear 
setting: Theorems 3.5 and Proposition 3. 

i? 8,2 regularity for the well-posedness in (1.1) is not sharp, but holds under very 
general assumptions on the coefficients and no smallness on data. While preparing 
this work for the publication, I have learned of a parametrix-based approach to the 
well-posedness of Quasi-linear Schrodinger for the small data in [18] in the lower 
regularity Sobolcv spaces than in the [14]. Adapting this approach to (1.1) may 
allow to lower regularity for small data. 
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The argument of the proof of the Theorem 1.1 also shows that (1.1) has a Local 
Smoothing effect, that is in addition to Xj, the solution is in u g LjH^Hx}^ dx). 

A simple transformation x — > —x in the equation (1.1) shows that the dispersive 
property (NL1) can taken with an opposite sign. That is, if in addition to (NL2)- 
(NL3), for every M > there exists a constant A m > 

»,,(.'•.(). < -A M |£| 2 

whenever (x,z) elx [— M, M] 3n , then Theorem 1.1 holds. 

We note that the continuous dependence in the Theorem 1.1 is the best we can 
hope for as (1.1) is quasilinear. 

Finally, the arguments in the proof of Theorem 1.1 can be used to extend the 
persistence of regularity to H s ' k rather than H s ' 2 . 

The rest of the paper is organized as follows. In the section 2 we prove the main 
linear estimate (1.3). In the section 3 we prove the well-posedncss of (1.6) for a 
time independent of the regularization e. Finally, in the section 4 we construct the 
solution of (1.1) and prove continuous dependence. 

Notation. When estimating with multiplicative constants, we often write A < x ,y 
B, to mean A < C(x, y)B, where the constant C(x, y) may change from line to line. 
As dependence on the integers (8,2) and the constant in the Sobolev embedding 
if 1 (R) c — > L°°(R) occurs frequently we do not explicitly state dependence on them. 



2. Linear estimate. In this section we prove the estimate (1.3) for the linear equa- 
tion (1.2) and show the necessity of several coefficient assumption. Main estimates 
for the well-posedness in sections 3.3 and 4 are later reduced to it. 

Let T' > and suppose (1.2) is valid on [0,T']. We assume that there exist 
constants A > 0, C\ > Co > and S > |, such that the coefficients a, b, c, d of the 
equation (1.2) satisfy the following conditions on R x7 = R x [0, T']: 

(LI) Dispersive property. Assume that the top coefficient is symmetric, that is 
a,ij(x,t) = ctji(x,t). Moreover, it is uniformly positive definite with 

<%(*, 0)6& > A|C| 2 

uniformly in x in R and £ € R". 
(L2) Regularity, a g C\B\, b G C}B 2 X , c g C\B\ and d g C}L™ with the space B k 
of bounded C k functions have the following bounds: 

||o(ar,0)|| B a + ||6(x > 0)|| Ba + ||c(x,0)|| B i + ||d(x,0)|| £ o. < C 

IhOM)!!^ + \\ b ( x >t)\\ciB2 + \\ c {x,t)\\ c i B i + \\d(x,t)\\ c i L ^ < Ci 
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(L3) Asymptotic flatness and decay. 

\\(x) 2S ' d x a{x,0)\\L- + \\(x) 2S b(x,0)\\L°° < Co 
\\(x) 25 d t d x a(x,T)\\ L ~ i+ \\(x) 2S d t b(x,t)\\ L ~ x 1 < d 
\\(x) s d x bHx,0)\\ L ~ + \\{x) s c\x,0)\\ L ~ < C 
\\{x) S d t d x b\x,t)\\L ¥x f \\(x) S d t J(x,t)\\ L ~ xi < & 
where and are the antisymmetric parts of b and c respectively, defined 

Theorem 2.1. There exist constants A = A(C , A, 5) and T = T(C , Ci, A, 5) < 1, 
such that if u(x,t) is a solution of (1.2) on I = [0,T'] C [0,T] , then u satisfies 
(1.3), i.e. 

\\u\\ L?Ll < A(\\H (x)\\ L 2 + \\f(x,t)\\ L]Ll ) 

and (1.4). 

By a solution of (1.2) we mean a classical solution and hence by the Sobolev 
embedding C^H 5 n CjH 1 would suffice. 

Regularity of a-d, particularly with weights in (L3) determines the Sobolev ex- 
ponent of H 8 - 2 in Theorem 1.1, where for simplicity, we set 6 = 5' = 1. Proof of 
analogues of Theorem 2.1 with coefficients rougher than above leads to lowering 
regularity in Theorem 1.1. 

Note, that for the applications of Theorem 2.1, the constant C\ will depend on 
the solution of the non-linear problem, while Co will only depend on data. As we 
will use the constant A = A(X, Co) from Theorem 2.1 to control the size of the 
solution and in turn Ci, it is crucial for A not to depend on C\. 

It is not difficult to prove an H s version of (1.3) and (1.4), provided coefficients 
arc more regular than (L2), by differentiating (1.2), using the Theorem 2.1 and 
choosing T small to control lower order terms. However, for some estimates in the 
proof of Theorem 1.1 we will not have such control of the coefficients, and we omit 
the H s estimate. 

Remark 1. For e = 0, Theorem 2.1 is valid with / = [-T',0] or [-T \T'] and 
CjH 1 n C°H 4 regularity. 

Remark follows from Theorem 2.1 by a simple scaling of the equation. Indeed, 
sending (x,t) — > (—x, —t) preserves assumptions (L1)-(L3), while changing the sign 
of the time. 

Finally before proceeding with the proof of the Theorem 2.1 in the subsection 
2.3, we motivate the coefficient assumptions (L1)-(L3) by showing the necessity of 
symmetry of a and decay of b for the Theorem 2.1. 

2.1. Symmetry of the top order. Similar to the hyperbolic systems, symmetry 
of (1.1) is necessary for the well-posedness. Namely, we consider the following 
constant coefficient linear system that violates the symmetry in (NL1), but satisfies 
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all other assumptions: 



u(x,0) T = (uo,i(x) u . 2 (x)j 



Taking a Fourier transform in space this equation reduces to an ODE. for which 
the explicit solution is 

tf = (e^^iff) + e^ 3 *fio,2(0) e« 8 *u , 2 K)) 

We then take the data uq(£) t = ( (^)~ s_ ) an d a computation shows 
||uo||if» = while for any t ^ 0, = oo. Therefore, (2.1) is ill-posed in 

H s for any s and hence in any H 11,2 . 

2.2. Necessity of decay of the coefficient b. Here we show that a slightly weaker 
form of asymptotic flatness (L3) is necessary for (1.3) in the following special case. 
More precisely, we show that for a solution u of 

{d t u + dlu + b{x)d 2 x u = f 
]u{x,Q) = uo(x) 

the condition on the coefficient b 

sup / b(x + r ■ s) rds < oo (2.3) 

(x,t,r) J a 

is necessary for (1.3). 

We show necessity by a WKB method, similar to an argument of [19] for a 
Schrodinger equation. Let u = e*^ x, *'^v(x, t, £), where 4> = x^ 2 + ££ 6 . Then for an 
operator L = d t + d x + b(x)d 2 we compute 

e^+UfPv) = (d t v - 3^d x v - b(x)^v) + [3if 2 Sg« + d 3 x v + 2b(x)i£ 2 d x v + b(x)d£v] 
We set v to be a solution of the transport equation 

f d t v - 3£ 4 d x v - b(x)^v = 
|u(x,0) = v (x) 

For which we get the explicit solution by the method of characteristics 
v(x, t, = e/o b (x+^^ id ° Vo (x + 3£ 4 t) 

Wc further define 

/ = e-** [3i£ 2 d 2 v + d 3 x v + 2b(x)i£ 2 d x v + b(x)d 2 v] 
Thus u solves (2.2), ||m||l 2 = MIl 2 an d ||uo||l 2 = ||«o||l 2 - 

Now assume that (2.3) does not hold. Then there exists a xq, to and ro such 
that 



I b(x Q + r • s)r ds > 61og(3A) 
Jo 
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Moreover we can assume by rescaling to with ^ that ro = 1. Hence for x near xq 

\ b(x + s)ds > 6 log(2A) (2.4) 
Jo 



Now define ig = ^ and likewise for any t. 

v(x, t*. = e'o 5 ^+ 3 « 4s )^o(z + 3£ 4 i«) - Jo ^+ s ) ds Uo (x + t) (2.5) 

Now let vq G C^°(K) be a function compactly supported near — to with IjuolU 2 = 
1. Then 



^ / (2^)^^(1)1^ = {2AY 

Moreover, by (2.5) and if b in £> 3 (i.e. C 3 , bounded and with bounded derivatives) 
we show that v(x,t^,^), . . . d x v(x, , £) are uniformly bounded in L x for G = 

|4|,|4|1. Therefore, ||/|| L j £ , < fh.^dt = o(l), as £ -> oo. However, 
(1.3) and calculations above imply that 2A < \\u\\LfLl < ^(1 + which is a 

contradiction for £ large enough. 

2.3. Proof of Theorem 2.1. To proceed we first change the dependent variable 
(or gauge it) and then proceed with the energy method for the gauged problem. 

We use (x)~ 2S to define a multiplicative change of variables as follows: 

<j>(x) = —N [ X (x'y 2S dx' with N = 2 + 10nd ° (2.6) 

Because the integral above is convergent and \d"(x)^\ < C(a, I3){x)^~ a we get 
<\> G B°° , i.e. it is bounded and all of its derivatives are bounded. Moreover, 

d x 0(x) = -N(xy 2S <0 

We now define the gauged variable 

v{x,t) = e~^u{x,t), (2.7) 

Definition of <f> implies that and e~^" are also in B°°, and hence 

NU 2 = ||e _0 -u|| L 2 Rijvi \\u\\lI 



Likewise, by the product rule and duality, 

||(a;) fc zi||_H- 3 ~N,s,k \\(x) k v\\ H 



" ,fc \\{x) k (d x ) s v\\ L i (2.8) 



for s = ±1. 



Inverting (2.7) we write (1.2) as: 

dtie^v) + a(a;,t)flg(e* (x >tf) + b{x,t)dl{e^ x) v) + c{x,t)d x (e' l ' {x) v) 
+ d(x, t)e^ x) v = -sdt{e^ x) v) + f(x, t) 
Then v satisfies the following "gauged" system 

d t v + Lv = —ed x v + eR^v + g{x, t) 
v(x,0) = v (x) 



(2.9) 
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where 

g(x,t) = e-+V)f( x ,t) and v (x) = e'^U^x) 
and the operator L = a(x, t)d x + b(x, t)d 2 + c{x, t)d x + dl has the coefficients 
bij(x,t) = bij(x,t) +3aij(x,t)d x (f)(x) 

Cij(x,t) = Cij(x,t) + 2bij(x,t)d x (t>(x) + 3aij(x,t)([d x (f>} 2 + dl<j>) 

dij(x,t) = dij(x,t) + Cij(x,t)d x </>(x) + b i3 (x,t)([d x cj)] 2 + d x <j>) 

+ aij (x, t){[d x (t>} 3 + 3d^d x (j) + d 3 x (b) 

and i?3 = [— d x ,e^] is a third order differential with B°° coefficients with norms 
controlled by N and 8. 

We now proceed with the energy estimates. Taking a dot product of (2.9) by v 
and integrating in x we get 

J dtv ■ v dx + J Lv-vdx = —e J d x v-vdx + e J R^v-vdx-\- J g-vdx (2.10) 

Note that L*, the adjoint of the operator L, is L*v = —a T d x v+l.o.t., where afj = ctji 
and l.o.t. are terms having less than 3 derivatives of v. Hence integration by parts 
of (2.10) gives 

^d t \\v\\ 2 L 2 = J Bd x v-d x v + Jcd x v-v + J Dv ■ v 

6 (2.11) 
— e / d x v ■ v + e I R3V ■ v + g ■ v = Ij 



3=1 



where 



/'.. = Bij = -3d x dji + bij + bji + (>(),0(iij 

- Cji = Cij = d x (bij - bji) + (cij - Cji) + 2(bij - bji)d x <f> 

II , = + d x bji - d x Cji + dij + dji 

We now claim that the Theorem 2.1 reduces to the following proposition. 

Proposition 1. There exist A = A{Cq, A, 8) and T = T(C\, Co, A, 5), such that for 
< t < T' < T < 1, the solution of (2.9) satisfies: 

d t \\v{x,t)\\l i + J(xy 25 \d x v\ 2 dx<A\\v(x,t)\\ 2 Ll +2\ J g-vdx\ (2.12) 

Indeed, discarding the second term on the right hand side of (2.12) and the 
Grownwall inequality imply 

\\v{x,t)\\li <e At (\Wo\\h+ J\ J g-vdx\dt) (2.13) 

Thus by the Cauchy-Schwartz inequality we get for / = [0, T'\ 

< e At (\\v \\ L2 + \\g\\ L]L2 J 

which is (1.2) after we use the comparability of the norms of u and v. 
Integrating (2.12) and using 

g ■ vdx = / {x) S {d x )~ 1 g -(d x ) ((x)~ S v)dx 
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in (2.13) gives the (1.4) after the comparability of the norms. 

2.4. Coefficient estimates. Before proving the Proposition 1, we first rephrase 
(L1)-(L3) in the following form 

Lemma 2.2. There exists T = T(C$, A, C\) small enough, such that for < t < 



L'l: 



uniformly in x. 
L'2: 



||aOM)|| c o B 3 + \\b(x,t)\\ c o B 2 + \\c(x,t)\\ c o B1 + \\d{x,t)\\ c o Laa < 2C 

L'3: 

\\(x) 25 d x a(x,t)\\L~ + \\(x) 2S b(x,t)\\ L ~ <2C 
\\{x) s d x tf{x,t)\\ L ~ + \\{x) 5 J(x,t)\\ L ~ < 2C 

Proof of Lemma 2.2. By the Fundamental theorem of Calculus, (LI) and (L2) 



&&>A|£| 2 -nTC 1 fc| ! 



a,ij(x,0)+ d t a i:j (x,t')dt' 
Jo 

Which implies (L'l) for T small enough. 

(L'2) and (L'3) are established similarly. □ 

Lemma 2.3. With <t <T as in Lemma 2.2 

(x) 2S By(x,t)6^ < -2ICI 2 (2.14) 

Moreover, there exists a constant A = A(Cq, \,S) > 0, such that 

\\(x) S C(x,t)\\ L ~ < VI (2.15) 
\\D(x,t)\\ L c* < A (2.16) 

Proof of Lemma 2.3. By the definition of B in (2.11), (L'l) and (L'3): 

(x) 2S Bijix^)^ < 6(x) 2S d x (j)(x)a ij (x,t)^ j 



+ n(x) 2 sup(3|9 a; a ii (x,t)| + \b i:j {x,t) + bji(x, 

ij 

< (-3N\ + l0nC o )\Z\ 2 



Which proves (2.14) by the choice of N in (2.6). 

Definition of C in (2.11), (L'3) and (2.6) imply (2.15). (2.16) is done similarly. □ 

Proof of Proposition 1. We estimate (2.11) term by term in reverse order: 
Trivially, 

h<\ v- gdx\ 
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For J5, by the Calculus of Pseudo-Differential Operators \d x ) 2 R3 is a ^DO of order 
I and hence maps to L 2 by the boundedness of Pseudodifferential operators. 
Thus 

h = £ / i? 3 « - (a.) 1 vdx < N , S e|MI^| 
Hence interpolation of ifi between L 2 and if 2 , Cauchy inequality and (1.5) give 

h < e\\v\\ 2 H i + A\\v\\ 2 L2 < e J \d 2 x v\ 2 dx + A\\v\\ 2 L2 
While for I4, integrating by parts gives 

h = — e I d\v ■ vdx — — e I \d 2 v\ 2 dx 



We use (2.16) for 1= 



h= I Dv-v<\\D{t)\\ L M\l 2 <A\\v\\\ 2 



By Cauchy- Schwartz and (2.15) we estimate 

(x) (x) 



/^|.(^C r «<||M| Us || (:c) ^T (t)||ioo||u1U2 

% % + A\\vf L , 
x) 



where Cf } = C j% . 



<\\-^\\ll+A\\v 

Finally, by (2.14) 

f, ,25 C^V , /" \d x v\ 2 , 

Summing together the estimates of 1\ - 1% we complete the proof of the Proposition 
1 and hence the Theorem 2.1 □ 



3. Wellposedness of the regularized problem. By Duhamel principle, solving 
(1.6) in a sufficiently regular Sobolev space is equivalent to a fixed point of the 
operator 

ft 

Yu £ ee e - Et9 *u + / e- £(t - t ' ) ^(-/V fl *u e + f)(t')dt' (3.1) 
Jo 

where the parabolic operator e~ * is defined as a Fourier multiplier: 

Well-posedness of (1.6) for a short time dependent on the parameter e is very 
similar to [14], where the same regularization is used for a quasi-linear Schrodingcr 
equation. However, as (1.1) is of higher order than Schrodinger , we provide the 
proof in the Proposition 2. 
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3.1. Preliminary estimates. 

Lemma 3.1. For any positive integer s, < T < 1 and t S I = [0, T] the following 
estimates hold 

II e- etd iu 1| c o H s,2 < C(s) \\u 1| Hs , 2 (3.2a) 

e 3 



e~ 



^- 4 '^ /(t')dt'||c ? H^ < T + d-\ ||/|| c o ff .-3, a (3.2b) 



By changing the Y norm in (1.8) by a multiplicative constant, for the rest of the 
paper we treat the constant C(8) as 1 from the Lemma 3.1 for simplicity. 

Proof. By Plancherel and boundedness of e~ a for a > 0: 



" et9 ^o|| 2 ff3 = / (0 2s e- 2et « 4 K(0! 2 ^ < No||^ (3.3) 



To proceed with weights, note that multiplying the PDE for v — e et9 *v,o by (a;) 
and commuting derivatives with weights, we get 

f d t (x) 2 v = ~edi ((x) 2 v) + E 3 ((x) v) 
\(x) 2 v(x,0) =(x) 2 u (x) 

Where £3 = [ed x ,(x) 2 } (x) 1 is a differential operator of order 3 with B°° coefficients. 
Hence by the Duhamcl formula 

rt 

(x) 2 e~ £t9 'u = e- £t9 i(x) 2 u a + / e^ 4 "*'^ E 3 {x) vdt' (3.4) 

Jo 

Therefore, by (3.3) and Minkowski inequality 

\\(x) 2 e- Et9 iu \\ H s < \\(x) 2 u \\ H , + [ We-^-^E^vWnsdt 1 



< \\(x) 2 u Q \\ H , +t\\E 3 (x) v\\ H s 

By the boundedness of (Pscudo)Diffcrential operators, \\E 3 (x) v\\h s \\(x) u||#s+s. 
By Pseudo-Differential calculus (or commuting (x) with derivatives by hand) and 
Cauchy-Schwarz 



IK*) v\\ H , +3 = j (x)(d x ) s (x)v-(x)- L (d x y +o (x)(d x yvdx < s \\{xy f\\i 3 \\f\\i 3+e 

(3.5) 

Thus by (1.5) 



\e- std ^u \\ H .,2 = \\v\\ H .,, < s \\u \\ H - 



which is (3.2a). 



Using elementary Calculus estimate a i e a < - for a > and the explicit 
formula for the semigroup e~ e x , t > we get 

* e-'W* ftf)dtf\\c° H . <s (t+ \\f\\c°H— 

Using this estimate instead of (3.3) with weights finishes the proof. □ 
We also need the following estimate for N(u) = Nu(u) from (1.7): 
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Lemma 3.2 (Moser estimate). Let u, v £ CjH s ' 2 for s = 8 or 9 be given with 

max{||t?|| c .o H 8,2, H|c°.ws>=} < M 

Then 

\\N{u)\\ c o H3 - 3 ,, < C s+3M (l + M s+3 )\\u\\ c » H sa (3.6a) 
\\N(u) - N(v)\\ c o Hs - s , 2 < C s+3 ,a/(1 + AI s+3 )\\u - v\\ c o Ha . 2 (3.6b) 

Proof. The proposition follows by the elementary calculus and the Sobolev embed- 
ding. The constants in the s = 9 case depend on M , because it is impossible to get 
terms like <9| +5 u ■ d^. +6 u by differentiating N(u) s + 3 times. □ 

3.2. Short time well-posedness. We now set up the following notation: 

AI = (16A + 2)R and T £ = min (— £ \\ (3.7) 



[^(l + M 11 )]^ 2 

with R from (1.8) and A from the Theorem 3.5. 

We do the contraction argument in the following closed subset of CjH 8 ' 2 for 
I = [0,T]. 

Xf = {u£C° I H S 2 : u(x,0) = fib(ar), \\u\\ c o H ^ < M, \\d t u\\c«H^ < C{M)} 

(3.8) 

with ||u — v\\ X " = \\u — v\ \c°h 8 ' 2 an< i 

C(M)= ^ sup \\f-N ri (u)\\ c o Hi , 2 +M 

11^11^0^8,2 < M 

Where by (1-8), (3.6a) and (3.7) C(M) is locally bounded and increasing in M. 

Proposition 2. T from (3.1) is a contraction map on Xf^ with I e = [0,T e ] and 
hence (1.6) has a unique solution u e in Xf 1 . 

We note, that (3.1) is locally well-posed in C t °i? 4 using the same proof. However, 
for the arguments in section 3.3 and beyond we need u e to be in Xf^ . 

Proof. Let u £ Xf 1 . Then by Minkowski inequality and (3.2), followed by (1.8), 
(3.6) and (3.7) we get: 

W)|| C[ oTe] tf 8 , 2 < IKIUb.s + \\f\\L^ Te] H»,z + (Te + ^)\\Nu\\ C ° ot] H^ 

M !~T~ 
<^r + (Ts+ y-j)C n , M (l + M )M < AI 
2 V £ fi 

Differentiating (3.1) in t and using ||ru(t)||c7° ^ ^ H s,2 < AI shows that ||<9irw|| c o^4.2 < 
C{M) and hence T(X^) C Xf^ . 

Likewise using Lemmas 3.1 and 3.2 for the difference gives the contraction prop- 
erty. □ 

Corollary 1. As we proved Proposition 2 by the contraction mapping argument, 
we automatically get continuous (and Lipschitz) dependence on data. That is, the 
flow map (mo, /) — > u is continuous in Y — > Xf*. 
Moreover, (3.1) has a persistence of regularity property. 
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For the persistence of regularity, let the data (uo,f) in addition to satisfying 
(1.8) satisfy for some s > 8 

IMIffO + H/IU^ ^ff*,2 + \\f\\c°_ 1A] H°-*-* < oo 

then the solution of (3.1) u £ £ XM satisfies ||w e ||r' H^-^nr 1 »■«— m < °°- To sec 

1 [0,T'] [O.T'J 

this, redo the boundedncss argument in the proof of the Proposition 2 for H 9,2 parti- 
tioning [0, T] into identical intervals of length Tg = T(e, M, C\ 2 ,m) for using (3.6a). 
For higher norms proceed inductively redoing (3.6a) to estimate H^Hco T] ff s - 2 usm g 

intervals of length f s = f(e, M s ,C s+3 A ~ u ) for M s = ||w e || c o H s-i,2. 

3.3. Adapting a linear estimate. We now aim to adapt the Theorem 2.1, in 
order to estimate u e uniformly in s. As the Theorem 2.1 applies to linear equations, 
we fix the coefficients (1.1) at an arbitrary v £ X^, to get 

d t u + N.;ju = —ed x u + f 
u(x, 0) = Uo(x) 

with N$ from (1.7). We will show in Lemma 3.3 that the coefficient bounds (Ll)- 
(L3) for this linear equation do not depend on the particular v, but only on the 
(NL1)-(NL3) bounds, data for t = and the bounds M, C{M) for t ^ 0. 

In particular an application of Theorem 2.1 for u = v = u E would imply for 
T' <T 

m\ L f Ll < A{\\u4 L , + \\f\\ L i Ll ) (3.9) 

We then pursue the same strategy for d x u £ and (x) 2 d x ~ 6 u £ for s = 6, . . . 14. 
Namely, wc differentiate (1.1) s times and account for quasi-linear interactions. 

' d t d s x tf + a (x, t, (d°u s ) a < 2 ) ■ d s + 3 u £ + b s {x, t, (d^) a < 3 ) ■ 8 s x + 2 u £ 
+c s (x, t, (d£u%< 4 ) ■ d^tf + d s (x, t, (d^) a < 5 ) d s x tf = -ed s + 4 u e , , 
+^/(i,i)+f s (i,i,(9 I "« £ )„< s _i)onlx/ l ' ' 

u £ (x,0) = u {x) on Rx{0} 

with a the same as in (1.1), 

6 s (x,t, (<9"u £ ) Q < 3 ) = sd x [a(x,t, [d x u e ) a < 2 ] + dg2a(x,t, (d x u e ) a < 2 ) ■ dlv 6 

+b(x,t, {di*tf) a < 2 ) + d^b(x,t, (d^u e ) a < 2 ) ■ d 2 u E 

where z = (?*) Q <2 = (d x u e ) a < 2 . Similarly, c s explicitly depend on a number s, 
(9"w e ) a <4 and up to 2 derivatives of a, &; and d s depend on s, (d x u 6 ) a <5 and up 
to 3 derivatives of a-c respectively. The reason for these formulas, is that differen- 
tiation is " linearizing" , so we only can get ad x +3 u e by applying all derivatives on 
u e and there are only few nonlinear ways to get terms higher than order s. 

Once we establish (L1)-(L3) coefficient estimates for (3.10) that are uniform when 
evaluated at v <G X^f in Lemma 3.3, then for the solution u e satisfying 

iT e CV ff 8+5 < 2 n c}H i+5 > 2 (3.ii) 

We would then apply Theorem 2.1 to (3.10) to get 

||a^|| c o ;L2 < A (||«U= + WlUj, i2 ) + AT'\\F s \\ LJiL , (3.12) 
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Provided T' < T. Note that, (3.11) is needed for to ensure that (3.10) is valid 
classically for 6 < s < 14. 

Finally, we multiply (3.10) by (a;) 2 to rewrite it as 

d t ((x) 2 d&) + adl{{xf d^) + b s d 2 ({x} 2 d° x tf) + c/d x ((x) 2 £%F) 
+d s (x) 2 d x tf = ~edi((x) 2 d x u*) + (x) 2 F°(x, t,tf,..., d°- l tf) 
+(x) 2 dsf(x,t)+EI(x)d°tf 
(x) d x u £ (x, 0) = (x) 2 d x u (x) 

where a, b s , c s , d s , F s are identical to (3.10) by construction and E 3 = [ed x + ad x + 
b s d x + c s d x ,(x) 2 ](x}~ 1 is an order 3 differential operator with L°° coefficients. Hence 
as long as the coefficient estimates are established and (3.11) is valid, Theorem 2.1 
implies 

\\(x) 2 d s X\\ c o /L2 <A(\\(x) 2 d s x u \\ L i + \\(x) 2 d s J\\ L]iL2 ) 
+ AT' (\\(x) 2 F s \\ LlL2 + \\E 3 (x) d s x ff\\ LlL2 " 



(3.14) 



provided T < T. 

Therefore, establishing a uniform estimate for H^H^m reduces to the following: 

• Find uniform bounds on the coefficients to (1.1) and (3.10) to use the Theorem 
2.1 (note the remark after Theorem 2.1). 

• Control the terms not involving data (F s and £"|) by choosing T small enough. 

• Remove the extra assumption (3.11). 

3.4. Uniform estimate. 

Lemma 3.3. Let u G X^r and consider the equation (3.10) for < s < 14 with 
coefficients evaluated at u, i.e. 

d t d s x u e + a (x, t, (d%u) a < 2 ) ■ s x +3 u £ + b s (x, t, {d%u) a < 3 ) ■ d s x +2 u e 

+c s (x, t, (d^u) a < A ) ■ d s x +1 u e + d s (x, t, {d*u) a < 5 ) d s x u e = -ed s x +A u e + f 

where b° = b, etc. Then the coefficients of this equation satisfy assumptions (Ll)- 
(L3) from the Theorem 2.1 with parameters A = X(R), Co = Co(s, R, C 3 ,r, Co) and 
Ci = Ci(s,M,C 3 ,m,C ,C(M)) with A(M), Cj M and C from (NLl)-(NLS), and 
C(M) from (3.8).' 

Proof. To simplify notation, let v = (u,d x u,d x u). Then ||w||c° i/»- 2 ~8,2 |Hlc°ifs,2 
and likewise for vq = v(x, 0) and dtv. Note that (1.8) becomes ||?7o||_ff 8 . 2 ^ R- 

By (NL1) a,ij(x,t,v) = aji(x,tv) and a(x,0,vo) > X(R) as by the Sobolev em- 
bedding HvoWIU 00 ^5 II ^o 0*0 Hi? 3 < R. Hence a(x,t, v) satishes (LI). 

To verify (L2) we use (NL2), Sobolev embedding and ||u ||^fs,2 < R: 
3 3 
J2\\diia(x,0Mx)}h~ < £ \\d2 d§a(v )d?v ---d2 im toh~ 

j=0 7oH \~H0\=O 

< C 3 ,h(1 + R 3 ) 
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Proceeding analogously and using (3.8), we get 

3 

Y,\\dtdla(x,t,v)\\ L ~ < C 3 ,m(1 + A/ 3 )(1 + C(M)) 

3=0 

and likewise for b s -d s . 

For (L3), we compute 

||(a;) 2 d x [o(x,0,iTo)]||i« < ll(^> 2 (d x a(x,0,v o ) + d s a(x,0,v o ) ■ d x v )\\ L ^ 
Using the Fundamental Theorem of Calculus 

d x a(x, 0, vq) — a(x, 0, 0) + / d x dga(x, t, 6vo)d9 ■ Vo 

Jo 

we get using (NL2), (NL3) and Sobolev embedding 

\\(x) 2 d x [a(x,0,v )]\\ L ~ < \\(x) 2 d x a(x,0,0)\\^ 

+ sup \\d^d^a(x, 0, g)\\ L oo \\(x) 2 d? 3 v \\ L - < Co + C 2M R 
|^l<ll«b||i«>;oi<l 

Similarly, 

\\(x) 2 d t d x [a(x, t, v{x, t)]||i- , < Co + 2C 2 , A/ (1 + M)(l + C(M)) 

Verifying (L3) for b s -d s is similar and involves Taylor expansion of coefficients at 
z — 0, chain rule and Sobolev embedding. □ 

Lemma 3.4. Let F s be from (3.10) and £J from (3.13) /or u e G X£{. Tften /or 
< s < 8 

||F s+6 (x,i, (c£<r) a<s+6 )|| if i2 < C(s + 6)C S+6>M (1 + M a+6 )||tr|| ifff . + « 
||(x) 2 F s (a;^,^^) Q<s )|U ri2 <C{s)C s , M {l + M s )\\(xftf\\ LTH s 
\\EI(x)d^\\ LTL 2 <C(s)C sM (\ + M%\\{x) 2 d^\\ LTL , + \\d^\\ LTH6 ) 

Proof. The first two estimates follow immediately from the construction of F s and 
Sobolev embedding. 

For the second, by construction E3 is a differential operator with coefficients 
bounded by C(s)C S: m (1 + M 4 ) and hence 

\\E s 3 (x) d^\\ LTL 2 < C(s)C s , M (l + M 4 )||(x) d s x tf\\ L?H s 

< C(s)C s , M (l + M 4 ) (\\(x) 2 d s x tf\\ L?L * + HfigtriU-H.) 

Where the last line follows by interpolation in (3.5). □ 

Theorem 3.5. There exists a constant A obtained using the Theorem 2.1 and 
depending on data, more precisely A = A{R, Xr, Co, Ci^r), and 
T = T(M,X R ,Co,C 14 , t M,C(M)) such that for any I C [0,T] a solution u £ G Xf 1 
of (1.6) on I satisfies 

\\u e \\ c o H s,2 < 8AR (3.15) 
Moreover, we can extend u e to be in X^ T i and solve (1.6) on [0,T]. 
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Proof. Note, that the Proposition 2 holds on [T',T' + T s ] with data (u £ (T')J) at 
t = T', as long as \\u e (T'), f\\ Y < 4f, which will hold for any T £ I as long as (3.15) 
is valid. When this is the case, u £ can be extended to solve (1.6) on [0, T" + T e ] in 
X^ T , +T j and arguing inductively to extend u £ to / = [0, T}. Thus it suffices to 
prove (3.15). 

We first complete the proof under the assumption that (3.11) is valid. 

By the Lemma 3.3 estimates (3.9), (3.12) and (3.14) are valid. Adding these 
inequalities and using equivalence of norms (1.5) we get 

IKIIcoh^ < A(\\u \\ H s. 2 + \\f\\ L]H s,2) + AT\\F 14 \\ LTL2 
+AT (j2W( x fF S hrv + £PS<*> WUfi» J 

\s=l s=0 / 

Using Lemma 3.4 we get 

\\U £ \\C°HS.1 < A(\\U \\ H S,2 + \\f\\ L l H S,2) + ATC 14 ,m(1 + M 1S )\\U 6 \\ C H S,2 

Which after choosing T small enough proves (3.15). 

To remove assumption (3.11) we regularize the data and note that (3.15) is 
uniform for every ||(uo, /)||y < R- More precisely, we convolve (uq, /) with a mol- 
lifier Xm{x) = mximx), where \ S Cq°(R) and / x(x)dx = 1 to get (wo,m,/m) € 
5? x C® ^S" and lim m _j. 00 (ito, m , f m ) = (uq, f) in Y. Therefore, for m large enough 
||(wo,m, /m)||y < R and by the Corollary 1 the corresponding solution u £ n £ Xf 1 n 
CjH 8+5,2 n C}H A+5 ' 2 and converges to u £ in JT^ whenever both are defined. 

Applying Theorem 3.5 to u £ m , for which (3.11) applies, allows us to recover (3.15) 
for u in the limit as m — > oo, as the constants A and T don't depend on m. □ 

4. Removing regularization. We first construct solution u of (1.1) in a topology 
weaker than C?iJ 8 ' 2 n C)tf 4 < 2 . 

Proposition 3. There exists aT > and / = [0,T] and a sequence of solutions of 
(1.6) w e „ m A 7 M /or e„ ->• 0, sucft that u 6n -> u in C a H 7 ' 2 n C}H 4 ' 2 . 

Proof. Take a difference of solutions u e and u £ £ Xf 1 of (1.6) for < e' < e < 1: 

<9 t (n e - u £ ') + A^u e - N r yu £ ' = -sd*(u e ~ u £ ') + (e - e')d^' 
u £ ~u £ '{x,0) = 

Then we can write N^e u £ — N~ r j u £ = L~ c (u £ — u £ ) where by the Fundamental 
Theorem of Calculus L- E has the form 

= a(x,t, (dlu £ ),< 2 )d 3 x + ~b(x,t, (d^),< 3 , {8^)^)31 
+ c(x,t, (d x u £ ) 3 < 3 , (diu £ ')j< 3 )d x + d(x,t, (d x u £ )j< 3l (^if')^)/ 
with coefficients 

hj = hj(x,t, u £ ,d x u El d 2 x u e ) + d x u E > :k / d z 2a lk (x,t 7 (d x ((1 - 9)u e + 9u e >) a<2 ))d0 

Jo 3 
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and similarly for c, d. 

By an argument identical to the Lemma 3.3, for u, v £ Xf' 1 , L^^ satisfies (Ll)- 
(L3) with bounds Co, Ci dependent on the same parameters as in Lemma 3.3. 
Hence for an appropriate T Theorem 2.1 applied to (4.1) gives 

\\u e -u e '\\c°H^ < A\\{e-e')dtu e '\\ L]Ll <ATM(e-e') (4.2) 

We then use that \\u £ — u £ \\ c o H a,2 < 2M and for < 1 < 1 interpolate H 13 between 
L 2 and H s+6 C H s ' 2 to get ' 

\\u e -u e '\\ LT Hii < \\u s -u s '\\^ L 2\\u £ -u £ '\§ rH ,o < [(e-e')AT\&(2M)® 

(4.3) 

Hence by completeness of CjH^. 3 we can take a sequence e n — > to construct 
u = lim^^o u 6 ■ 

Note, that LfH li = (i}i/~ 14 )*, and hence the closed unit ball in it is weak-* 
compact by the Banach-Alaoglu theorem. As a consequence, the sequence u £n G 
Xf 1 C LfH^' 2 has a weak limit u £n LfH^ up to a subsequence. Hence, for 

a.e. x, there is a further subsequence, such that u £n (x) — > v(x). But by uniqueness 
of limits v — u. A similar argument for LJ°i? 8 ((a;) 4 dx) gives 

u £ » — ^ u G LfH 8 - 2 (4.4) 

Using Fatou lemma gives ||{t|| LO o ff 8.2 < M. 

For strong convergence in Cji? 8-1,2 we multiply (4.1) by (a;) 2 and rewrite it as 

d t (x) 2 (u £ -u £ ') + L^ ^, (xf [u £ -u £ ') = -ed 4 x (x) 2 (u £ - u £ ) 
+ (e - e')(xf d±u £ ' + E 3 (x) (u £ - u £ ') 
(x) 2 {u £ -u £ '){ Xl 0) = 

with i?3 = [— ed* — L- e ge ' ,(x) 2 ] is an order 3 differential operator with coefficients 
bounded by 4 + Ci,m{\ + M). Applying Theorem 2.1 all terms are treated as in 
the proof of (4.3), except we use (3.5) to interpolate 

\\E z (x) (u £ -u £ ')\\ LTLl < M , Cl , M \\u £ - iT \\l TH6 \\(xf (u £ - u £ ' )\\ LTL 2 ^0 

as e, e — > 0. Interpolating 7J 8_1 between L 2 and H 8 implies (x) 2 u e — > (x) 2 u in 

Finally, using (4.1) and Lemma 3.2 implies dtu £n —> dtU G CjH 4 ' 2 . □ 

Moreover, for uniqueness of the limit we use (4.3) for u £ = u and u £ = it' for 
e = e' = 0. The regularity in the Proposition 3 is enough for the Lemma 3.3 to be 
valid. Moreover, as Theorem 2.1 is valid for [— T, T], when e — we get uniqueness 
for [— T, T]. 

We now aim to recover the loss of the derivative in Proposition 3 following the 
regularization method of Bona-Smith [2]. That is we regularize the data with a 
parameter k and then apply the Theorem 2.1 for the difference at the top level of 
regularity keeping a careful track of n. 



20 



TIMUR-AKHUNOV 



0(KI) = 



4.1. Data regularization. Let < < 1 be a radial smooth bump 

1, if |£| < 1 

0, if iei > 2 

Define for all < k < 1, wo. K = ■ 0(k|£|). 

Lemma 4.1. Let K. <e H 1 * be a compact set. Then Vk > 0, and any uq £ fC, 
uq.k £ & satisfies 

\\uo,k\\h 14 +3 = 0(k~ j ) for all j > (4.5a) 
|juo : K - "o||l 2 = o(k 14 ) and ||u ,« - wo||h 14 = o(l) (4.5b) 
mi/i i/ie convergence rate dependent on K.. 
Proof. Let j > and < k < 1 

IKJSri** = / (i + l£l 2 ) 14 M£)l 2 (i + itfymtlM < (-) 2j h \\ 2 H ^ 
J\i\<- K 

(4.6) 

Which proves (4.5a). 

For (4.5b) is suffices to show the first estimate, with the second done identically. 
By the Fundamental Theorem of Calculus: 

/ \ 2 



Krt-UoHia = /(<M^I)-l) 2 M£)| 2 d£< / Sup \4>'(V)\ K 2 |^0| 2 ^ 

where we used 0(0) = 1 and defined B K \ ( \(0) = [-k|£|, k|£|]. As 9^(0) = for all 
j > 0, we can continue with the Taylor expansion of |<^^(t7)| 14 times and then use 
dg(f) = on B±(0) for j > to conclude 



uo\\h < K 28 \\diu\\u [ \e 4 uoM 



"0,, 

'i?i>- 

The o(l) rate comes by the Lebesgue Dominated Convergence that is uniform for 
uq in a compact if. □ 

Lemma 4.2. Lemma 4-1 is va/id weights. That is for K d if 8 ' 2 

\\{x) 2 u^ K \\ H &+j = 0{n~ r ) for all j > 

\\(x) 2 (u ,« - w )||l2 = o(k 8 ) and ||(:r) 2 (u 0)K - Uo)||#s = o(l) 

Proof. The Lemma 4.1 is valid for ((a;) 2 uo) K (with obvious modification of 14 to 8), 
so it suffices to show that v K = (x) uo, K — ((x) uq) k has norms = 0(k~ n ) 

for all s and N. 

Using (a;) 2 — (x — y) 2 = 2(x — y) ■ y + \y\ 2 we get 

(2(x ~y)-y+ \y\ 2 ) u (x - y)-<p{~)dy 

K K 

Hence 

and as d^4>(0) = for j > 1 we can get arbitrarily high power of k by Taylor 
expanding d^cft. Thus remark follows, as xuq and uq € L? . □ 
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Remark 2. Lemmata 4.1 and 4.2 hold for JC <e L\H 8 - 2 or /C <e C^H 8 ' 2 . That is 
for each / G /C, there is / K G such that 

- = o{n li ) \\{xf (f K f)\\ L]L2 = o(k 8 ) 

- /llijfl.." = WhhlH^ = O(K-i) for all j > 
and similarly for K <s CjH A > 2 . 

For L 1 !! 8 ' 2 redo the Lemmas using 

lim /(/ (l + |^| 2 ) 14 |/| 2 dO^^ = 0, 

for o(l) convergence rate, which follows from Lebesgue Dominated convergence in 
t and £. 

Finally, for / G CjH 4 ' 2 we use the continuity of / to see that /(/) is a compact 
set in H 4 ' 2 . 

Putting together Lemmata 4.1 and 4.2 and Remark 2 we get 
Proposition 4. 

\\(u ,f)-(u 0!K ,f K )\\ Y = o(l) 

and (u , K ,f K ) G H 8+s > 2 x (L|_ 1 ^if^ 8 ' 2 n C[ 1 _ 1 ^H^ 3 - 2 ) with 0(k- s ) norms for 
s > 0. 

We now want to improve the result of the Proposition 3 for the solution u e K of (1.6) 
for data (io, K ,/ K ). By the Proposition 4, for k > small enough ||(uo,re, /k)||y < 
R, hence the Theorem 3.5 applies to the u £ K . Moreover, by Proposition 2, u e K G 
CjH 8+s ' 2 n C}H 4+S ' 2 for any s > 0. In fact, u% satisfies the following uniform 
bound 

Proposition 5. For interval I from Theorem 3.5 there exist a constant C indepen- 
dent of e and k, such that for j = 1,2, 3: 

\\K\\c°HS+J-2 < C(||uo, K ||/f8+j.2 + \\f K \\ L i H «+j-2) = 0(k^ 3 ) (4.7) 

Proof. Let I k = [kf, (k + l)f] for f to be fixed below and < kf < T - f . 

Then considering (1.6) with data (u e K (kf) , f(t + kfu Lemma 3.3 is valid on Ik 

and Lemma 3.4 can be refined to hold for < s < 8 + j. Namely, there's no way to 
get terms involving d], A u e ■ <9j. 4+:, £t £ by differentiating N^{u e ) 17 times. This gives 

ll^llc? H«+i:* < C{\\u e K (kf)\\ H s +] ,2 + \\f K \\ L i H »+j.z) 

J k 1 

+Cf C 14+j , M (1 + M 18 +->) || if || c o H s +j , 2 

'k 

with C = A(|| (u% , / k )||f) from Theorem 2.1. Choosing T small enough and induct- 
ing on k finishes the proof. □ 

Hence interpolating \\u% — u% \\c°h 14 between L 2 estimate in (4.2) and H 15 in 
(4.7) gives that u K G CjH 14 . Doing the same with weights, implies 

u K G C°H 8 - 2 (4.8) 
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4.2. Uniform convergence. 

Theorem 4.3. u £ K — >• u £ as k —> in UfB\ 4 uniformly in < e < 1. Moreover, 
this convergence is uniform for u £ coming from data in a compact set K <e Br(0), 
where Br(0) is a ball of radius R centered at in Y . 

Proof. Redoing (4.2) for u £ and u% we get 

\\u £ - ul\\ LT Ll < A(\\u - u , K \\ L 2 + ||/ - f K \\ L ]Ll) = o(k 14 ) (4.9) 

where the last identity is by (4.5b) and Remark 2. 

We then subtract (3.10) for u £ K from (3.10) for u £ with s = 14 and rewrite it as 

d t dl\tf - u%) + N u (u £ - u £ K ) = -edl\u £ - u%) + dl\f- f K ) + F 14 - 
+ (N u .k - N u )u% 

dl 4 (u £ - u%)(x, 0) = dl 4 (u {x) - u , K {x)) 

where N\± is the operator from the Lemma 3.3 for s = 14 with coefficients evaluated 
at u £ and F 14 = F 14 (x, t, u £ , . . . , d 43 u £ ), and likewise iVu jK and F^ 4 are evaluated 
at u £ K respectively. Therefore, assuming (3.11), we apply the Theorem 2.1: 

WdlHu* - uDhfLi < A(\\dl\u - u , K )\\ L% + \\dl 4 (f- f K )\\ L]Ll ) 

+ A\\F li - F^\\ L]Ll + A\\(N U , K - N u )u £ K \\ L]Ll =h+h+h 

We estimate the right hand side term by term. By the Proposition 4, I\ = o(l) as 

hi — y 0. 

For li we write by the Fundamental Theorem of Calculus, 
f 1 r)F 14 

F 14 -P?=^2 -j— (x, t, (dP(9u £ + (1 - 0)O/*<i3) d9 ■ d«(u £ - u £ K ) 
Hence by Sobolev 

h < ATC l6 , M (l + M 14 )\\u e - t£|Uf a;« 
We then interpolate H 13 between L 2 and H 1A as in (4.3) using (4.9): 



For J 3 , we proceed as for I 2 using Fundamental Theorem of Calculus 




Hence by the Cauchy-Schwarz and Sobolev: 



\\(N 14>K - N u )u £ K \\ L]Ll < C 5)M (1 + M 3 )\\il £ - tfJ L?H »\\dl 7 ir K \\ L?Ii 
Thus by (4.9) and (4.7) 

h = o(k 17 )0(k- 3 ) = o{k 14 ) 

Combining the estimates for terms I1-I3, (4.9) and using equivalence of norms (1.5) 
gives u% u £ in CjH 14 provided (3.11) is valid. 

To finish the proof without (3.11), we proceed as in the proof (3.5), except 
when we regularize the data (uo im , fm) — > {uq, f) in Y, we use compactness of 
K = {(uo, m , f m )} U {(uo, /)} in Y to ensure that the convergence rate u £ m K — > u £ nl 
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as k — > is uniform in m. 

Finally, as the rate of convergence in k comes from Lemma 4.1 as longs as 
IK^o, j)\\y < R, it remains valid for data in a compact set K <g Br(0). □ 

Redoing the Proposition 4.3 with weights, we get u £ K —> u e in CjH 8 ' 2 . 

4.3. Time continuity and continuous dependence. By the Proposition 3 u e — 
u% u — u K in LfH 14 . Thus for fixed K, we apply Fatou's lemma to J d]. A {u e — 
u E K ) ■ <j>dx for a vector valued test function <fi £ Cq° (R) to get 

\\u-u K \\ L f>H li < liminf||M e - u E R \\ LrH u = o(l) (4-10) 

1 x 6— >0 1 x 

as k — V 0. Thus the Theorem 4.3 is also valid with e = 0. 

Combining (4.10) with (4.8), gives u € C^H 14 . Mor eover, 
u e -u = (u £ - u e K ) + (u K -u) + {ii s K - u K ) 

with the first two terms that go to as n — > in CjH 14 uniformly in e and the last 
term goes for each fixed k. Hence u e — » u in CjH 14 . 

Redoing the same argument with weights proves u e CjH 8 ' 2 

lim||u £ -u|| c o ff8 , 2 =0 (4.11) 

with the convergence rate dependent on the data profile, namely the compact set 
JC<sY. 

Proposition 6. Suppose (uo.m, fm) (uo, /) in 5^- TAen i/ie respective solutions 
of (1.1) u ro u in Xj = C^ 8 - 2 n C}H^ 2 . 

Proof. For m large enough the data have size less that R. For such m we write 

Wm - U = (U m - 4) + ~ u) + «„ ~ U £ ) 

The first two terms converge to in CjH s ' 2 as e — > uniformly in m by compact- 
ness of data, while the third goes to as m — > oo for each fixed e by the Corollary 1. 

For the C}H 4 ' 2 use (1.1) and Lemma 3.2. □ 

Finally, for Persistence of Regularity for (1.1), proceed inductively s = 15, 16, . . . 
as in the Proposition 5, applying Theorem 2.1 for (3.10), so that for each step terms 
from the Lemma 3.4 are controlled by the previous step. Then remove 3.11 by the 
Proposition 6. This completes the proof of the Main Theorem 1.1. 
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